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ABSTRACT 

Context. Tow important stellar-population parameters (age and metallicity) of the dominant stellar populations (DSPs) of galaxies are 
usually estimated by comparing the observed absorption line indices or colors to predictions of some simple stellar population models. 
However, some studies show that there is actually recent star formation in galaxies, including early type ones. This suggests that we 
may not be obtaining accurate the two stellar-population parameters for the DSPs of galaxies. This is obvious when we estimate the 
two parameters by colors, because the youngest populations dominate the light and make the fitted stellar populations younger and 
richer in metal. 

Aims. We plan to study how young populations (YSPs) in composite stellar populations (CSPs) affect the colors of star systems and 
to analyze how the stellar ages and metallicities derived from colors possibly deviate from those of the DSPs. 
Methods. The techniques of stellar population synthesis and Monte Carlo process are used in this analysis. 

Results. It is found that the age and mass fraction of a YSP affect colors of a mixed star system significantly, but the former is 
stronger. In addition, our results show that the stellar ages and metallicities derived directly from a pair of colors are about 2. 14 Gyr 
younger, while 0.0027 more metal rich on average than those of the DSPs of composite stellar systems. Some possible distributions 
of the differences between stellar-population parameters determined by colors and those of DSPs of CSPs are presented. The possible 
distributions of the differences between colors of CSPs and those of their DSPs are also shown. Stellar ages and metallicities measured 
by colors and line-strength indices are compared in the work, with a sample of 18 galaxies. Furthermore, the YSPs may affect the 
fundamental plane and Kormendy relation of early type galaxies. 
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1. Introduction 

To study the formation and evolution of galaxies, we usually 
need to estimate two stellar-population parameters (age and 
metallicity) of the dominant stellar populations (DSPs), which 
contribute most to the stellar mass of galaxies. The two pa- 
rameters of DSPs of galaxies without obvious star formation, 
e.g., early type galaxies, are usually estimated by comparing the 
observed values with the predictions of simple stellar popula- 
tions (SSP s) (see , e.g., Worthev 1994,: Vazdekis [19991 Bruzual 
& Ch ariot 12003] (hereafter BC03); Fioc & Rocca Volmerange 
[T997(PEGASE code); Zhang et al. '2005'). Most of these works 
are accomplished by measuring absorption line indices or spec- 
tra, e^Kuntschner (2000), Vazdekis et al. (1 19971 1. Gallazzi et 
al. (|2^ ). Li et al. (,200 ^, Zh ou et al. ( fT992] l. Kong & Cheng 
( I1998I I. and Kong et al. (I2003I I. A few works have also tried to 
obtain similar results by measuring colors, e.g., Dorman et al. 
( I2003l l, Yi ( 120031 ), Wu et al. (2005 ), James e t al. (2 Q06 ), Li et al. 
( I2007I I. Kaviraj et al. (2006 ), and Lee et al. (I2006I I. because col- 
ors are very useful for studying the stellar populations of distant 
galaxies. However, many studies show that there is recent star 
formation in galaxies of all types (Trager et al. 120001 Yi et al. 
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120051 Schawinski et al. I2006I I. which suggests that we have not 
obtained accurate ages and metallicities of the DSPs of galax- 
ies, via SSPs with instantaneous star burst (see, e.g., Serra & 
Trager l2007l Li et al. 120071 ). This is true when one uses colors 
to estimate the two parameters. If there are two populations in 
galaxies, we will usually get younger and more metal-rich stel- 
lar populations than the DSPs of galaxies. Thus it is difficult to 
get reliable estimates for the ages and metallicities of DSPs of 
galaxies from colors before we know how much the two val- 
ues derived from colors are possibly different from those of their 
DSPs. Unfortunately, there has been no detailed study of this yet. 

We intend to get some answers in this work. The main aims 
of the paper are to study how the parameters relating to the 
young population (YSP) of a composite stellar population (CSP) 
affect the colors of the system and to analyze the distributions of 
the deviations in stellar ages and metallicities determined by the 
colors of CSPs when comparing them to their DSPs. We also 
intend to study the possible distributions of the differences be- 
tween colors of CSPs and their DSPs. 

The structure of the paper is as follows. In Sect. 2, we briefly 
introduce the construction of CSPs. In Sect. 3, we analyze how 
the parameters of the YSP of a CSP affect the colors of the sys- 
tem. In Sect. 4, we study the distributions of the differences be- 
tween colors of CSPs and those of their DSPs. In Sect. 5, we give 
some possible distributions of differences between stellar ages 
and metallicities estimated via colors of CSPs and those of their 
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Table 1. Input parameters for modeling composite stellar popu- 
lations. 



parameter 



values 



[Z/H], 
U (Gyr) 
[Z/H]2 
h (Gyr) 

F2=M2/(Mi+M,) 



-2.30,-1.70, -0.70, -0.40, 0.0, 0.40 

0.05, 0.1, 0.51, 1.02, 1.61, 2, 2.5, 3, 20 

from [Z/H], to 0.40 

from 0.01 to fi 

0.005, 0.01, 0.015, 0.5 



DSPs. Then we try to study how we can correct for the devia- 
tions in the two stellar-population parameters, which result from 
the existence of YSPs in galaxies. The stellar ages and metal- 
licities estimated by different methods (the color-method and 
absorption-line index method, a color pair method and multi- 
color method) are also compared in the section. Finally, we give 
our discussions and conclusions in Sect. 6. 



2. Construction of composite stellar populations 

We construct CSPs using the SSPs of BC03, with a wide range 
in age (0.01 - 20 Gyr) and metallicity (0.0001 - 0.05). Each 
CSP is constructed of two SSPs: an old DSP and a YSP There 
are five input parameters for each CSP: metallicity and age of 
the old DSP (Zi and t\), metallicity and age of the YSP (Z2 and 
fa), and the mass fraction of the YSP in the total system (F2). 
For the convenience of comparing our results with others, we 
mainly use [Z/H] 1 =log(Zi/Z0) and [Z/H]2=log(Z2/Z0) to repre- 
sent the stellar metallicities of the DSPs and YSPs in this work. 
The range of F2 is from 0.5% to 50%. Note that f2 < f 1 , Z2 > Zi . 
In this way, we construct 1 894 200 CSPs and then calculate all 
their UBVRIJHK colors. The ranges of input parameters of the 
sample CSPs are listed in Table 1 . 



3. Effects of the YSP parameters on colors 

A young population usually makes the whole star system bluer 
than the DSP of the system and results in incorrect determina- 
tions of stellar age and metallicity. However, we do not know 
how the parameters of the YSP affect the colors of the whole 
system. 



3.1. Effect of the age of the young population 

We test the effect of the age of the YSP of a CSP, f2, on the 
colors of the CSP via mixing YSPs with variable ages into the 
same DSP. Then we investigate how the differences between the 
colors of CSPs and those of the DSP change with the age of 
the YSP. We take the same metallicity for the YSP and DSP of 
a CSP here. The results for three colors that have the potential 
to disentangle the stellar age-metallicity degeneracy (Li et al. 
I2OO7I 1. i.e., (B - V), (B - K), and (/ - H), are plotted in Figs. 
1, 2, and 3, respectively. The differences in the colors, which 
are represented by A(B - V), A(B - K), and A(/ - H), etc., are 
calculated by 



AC = Ccsp - C, 



DSP, 



(1) 



where AC is the difference, while Ccs? and Cdsp are the col- 
ors of the CSP and DSP, respectively. We take a few CSPs with 
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Fig. 1. Differences between {B - V) colors of a CSP and DSP 
versus YSP age {12). Solid, dashed, dot-dashed, and dotted lines 
con-espond to YSP mass fractions {F2) of 0.5%, 1 .0%, 5.0%, and 
10.0%, respectively. A(Z? - V) is calculated by Eq. (1). 



iog{*i/yr)=io-o (b) 

[Z/H]2=[Z/H],=0.0 




iog{i,/yr)=io.o (d) 

[Z/H]2=[Z/H],=0.40 




Fig. 2, Similar to Fig. 1, but for {B - K). 



log(f I /yr) of 10. 1761 (f I = 15 Gyr) or 10.0 (fi = 10 Gyr) and metal- 
Ucity ( [Z/H] 1 = [Z/H] 2) of 0.0 or 0.40 in this study. Four mass 
fractions (0.5%, 1.0%, 5.0%, 10.0%) are chosen for the YSPs in 
the tests. As we see, the results for different colors are similar: 
The differences between colors of a composite star system and 
those of the DSP of the system increase rapidly with decreasing 
age of the YSP when the age is less than about 5 Gyr. The results 
for some other colors and other CSPs, which are not shown in 
the paper, are found to be similar to those shown here. 

3.2. Effect of the mass fraction of the young population 

The effect of the mass fraction of the YSP, F2, on the colors of 
a CSP is studied by a method similar to what was used to study 
the effect of the age of the young population. There A{B - V), 
A(B - K), and A(/ - H) are shown to increase with F2. We plot 
the results in Figs. 4, 5, and 6, respectively. We take 15 Gyr 
(log(fi/yr)=10.1761) or 10 Gyr (log(fi/yr)=10.0) for the ages 
of the DSPs and take 0.51 Gyr (log(f2/yr)=8.7076) or 2 Gyr 
(log(f2/yr)=9.3010) for the YSPs here. In each test, we adopt 
metallicities ( [Z/H] 1 = [Z/H] 2) of -1.70, -0.70, 0.0, and 0.40 for 
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Fig. 3. Similar to Fig. 1, but for (/ - H). 



Fig. 5. Similar to Fig. 4, but for {B - K). 
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Fig. 4. Difference between (B - V) of a CSP and that of its DSP 
versus F^- Solid, dashed, dashed-dotted, and dotted lines show 
CSPs with metalHcities ( [Z/H] 2= [Z/H] i) -1.70, -0.70, 0.0, and 
0.40, respectively. A(B - V) is calculated by equation (1). 
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Fig. 6. Similar to Fig. 4, but for (/ - H). 
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4. Possible distributions of the differences between 
colors 



the CSPs, which are shown differently in the figures. The metal- 
licities of the YSP and DSP of a CSP are assumed to be the same. 
We see that the results for three colors are similar: (1) Compared 
to those of DSPs, colors change quickly with increasing mass 
fraction of the young population, especially for mass fractions 
lower than 10%. (2) The rate of change in the colors does not 
mainly depend on the difference between and fi, but on the 
age of the YSP 

3.3. Effects of the metallicity of the young population 

We also investigated the effects of the metallicity of the young 
population, i.e., [Z/H]2, on the colors of a CSP. We use a similar 
method to the one above. Some parameters, i.e., fi, [Z/H]i, f2, 
and F2 are fixed for CSPs, while [Z/H]2 are variable in each test. 
Note that we assume [Z/H]2 > [Z/H]i. The results for (B - V) 
are shown in Fig. 7. We see that [Z/H]2 or Z2 does not affect 
the colors of star systems as clearly as t2 and F2, but the effect 
is obvious for CSPs with metal-rich DSPs. The results for other 
colors are not shown here, since they are similar to those of (B - 
V). 



The difficulty of getting accurate stellar ages and metallicities 
for the DSPs of galaxies via colors is that we get the colors of 
CSPs rather than those of DSPs. In this section, we analyze the 
possible distributions of differences between colors of CSPs and 
those of the DSPs mainly via a Monte Carlo technique. 

4.1. The composite stellar population sample 

We construct a set of CSPs for the study via the following steps. 
First, we generate the ages and metallicities (fi and [Z/H]i) of 
DSPs in our sample CSPs according to two distributions de- 
rived from Gallazzi et al. ( I2005I I. who determined stellar ages 
and metallicities of 175 128 galaxies drawn from the Sloan 
Digital Sky Survey Data Release Two (SDSS-DR2) by com- 
paring D4000, HyS, HdA+HyA, [Mg2Fe], and [MgFe]' indices of 
galaxies to predictions of SSPs of the BC03. In particular, all 
galaxies have a ratio of signal-to-noise S/N > 20 and do not 
have current star formation. Stellar ages used in this work are 
corrected for the light travel time (Ho = 70 km s ' Mpc '). This 
will not obviously change our result, because the colors of CSPs 
are more sensitive to the age of the YSP, t2, rather than that of 
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Fig. 7. Difference between (B - V) of a CSP and that of its DSP 
versus Z2. Solid, dashed, dot-dashed, dotted, and dash-dot-dot- 
dot are for metallicities of DSPs ([Z/H]i) of -2.30, -1.70, -0.70, 
-0.40, and 0.0, respectively. 



the DSP, f 1 . However, the correction can increase the difference 
between fj and t2 slightly, and, correspondingly, the difference 
between colors of a CSP and its DSP. This makes the difference 
calculated in this work closer to the actual difference between 
colors of galaxies and their DSPs, as galaxies possibly contain 
more than two SSPs. Second, the metallicity of the YSP of a 
CSP is taken as the same as that of the DSP. Third, the ages of 
YSPs are given randomly by a Monte Carlo process, with a nat- 
ural assumption that the YSP of a CSP is younger than the DSP 
of the CSP (0 < f2 < fl). Fourth, the mass fraction (F2) of YSP 
is calculated according to the ages of the DSP and YSP by the 
equation 



h - 1\ 
Fo . exp [ ], 



(2) 



which assumes that the mass fraction of the young population 
declines exponentially with decreasing age of the YSP, accord- 
ing to some studies about the star formation histories of galaxies, 
e.g., Bruzual ( 1983) and Thomas et al. (|2005l l. The variables ti 
and t[ are the ages of the YSP and DSP of a CSP, respectively. 
The value of Fq is the maximum mass fraction, while t is a free 
parameter. For a standard case, we take 0.5 and 3.02 for Fq and 
T, respectively, which gives about 15% CSPs with F2 < 4%. This 
agrees with the result of Yi et al. (2005) that there are about 15% 
bright early type galaxies with recent (< 1 Gyr) star formation 
at a level more than 1% - 2% in mass compared to the total 
stellar mass. In this way, 30 000 CSPs with f 1 > 1 Gyr are con- 
structed. The distributions of the input parameters of these CSPs 
are shown in Fig. 8. 

4.2. Distributions of differences between colors of CSPs and 
fhose of tineir DSPs 

We investigate the distributions of the differences between colors 
of CSPs and those of their DSPs using the above sample. The 
results of eight colors, which are shown to be relatively sensitive 
to stellar age or metallicity (see Li et al. 2007), are plotted in 
Figs. 9 and 10. The differences are calculated by Eq. (1). We see 
that more CSPs exhibit bluer colors than their DSPs. 
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Fig. 8. Distributions of five input parameters of our sample 
CSPs. Panels (a), (b), (c), and (d) show the distributions of 
[Z/H]i and [Z/H]2, log(fi/yr), log(f2/yr), and F2, respectively. 
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Fig. 9. Distributions of differences between colors of CSPs and 
those of their DSPs. Panels (a), (b), (c), and (d) show the distri- 
butions for {B - V), (R - /), {B - K), and (/ - H), respectively. 
A(B - V), AiR - /), A(B - K), and A(/ - H) are calculated by Eq. 
(1). 

5. Possible corrections for the deviations in stellar 
ages and metallicities determined by colors 

5.1. Distributions of tliese deviations 

We investigate how the stellar ages and metallicities measured 
by colors deviate from those of the DSPs in the sample gener- 
ated in Sect. 4. Because the two stellar-population parameters of 
galaxies are usually estimated via two colors due to the limita- 
tion of the observational data, we study the results derived from 
different color pairs in the work. In detail, stellar ages and metal- 
Ucities are determined by [(B - V), (V - K)], [(B - V), (B - K)], 
and [(B - V), (/ - H)], respectively, using the BC03 SSPs. The 
three pairs of colors (hereafter color pairs) are selected by find- 
ing the best pairs for breaking the stellar age-metallicity degen- 
eracy (Worthey 119941 ) "by eye" via color-color grids of color 
pairs constructed via three metallicity-sensitive colors, (B - K), 
(I - H), (y - K), and thr ee age -sensitive colors, {B-V),{U - R), 
and {R - 1) (see Li et al. |20^ . Note that [{B - V), (V - K)] and 
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Table 2. Coefficients of Eqs. (3) and (4). 
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Fig. 10. Similar to Fig. 9, but for (U-R) (Panel a), (V-K) (Panel 
b), (/ - K) (Panel c), and (U - K) (Panel d). 

{{B - V), (B - K)] are equivalent from a mathematical view, as 
they relate to the same three magnitudes, but they have differ- 
ent abilities for disentangling the well-known degeneracy. For 
example, it seems impossible to disentangle the above degener- 
acy via [(B - K), (V - K)]. The reason is that both (B - K) and 
iV - K) are relatively sensitive to metallicity, but their relative 
metallicity sensitivities are different (Li et al. 2007 ). Thus it is 
necessary to compare the results derived from [{B - V), {B - K)] 
and [(B - y ), (y - ^:)] . Then the deviations ( AAge and AZ) of the 
best-fitted values (age, ff and metallicity, Zf ) from those of their 
DSPs (age, t\ and metallicity, Z\) are calculated (AAge = ff - 1\, 
NZ - Zf - Zi). The distributions of the deviations in stellar ages 
and metallicities are plotted in Fig. 11, in which the results re- 
ported by three different pairs of colors are shown, respectively. 
As we see, the results derived from fitting three different pairs 
of colors are similar. They have similar shapes and all of them 
peak near AAge = -0.8 Gyr and AZ - 0.0002, with averages - 
2.14 Gyr and 0.0027 in stellar age and metallicity, respectively. 
For the convenience of using these distributions in future works, 
we fit them by equations: 



P(AAge) 



f'o + ^Ai exp[ 

i=l 



(At-Xj)- 
-2W^ 



P(AZ).Po'.iA:exp[^^^], 

i=I i 



(3) 



(4) 



and the coefficients are listed in Table 2. Here P(AAge) and 
P(AZ) in ffie two equations are ffie percentages probabilities. The 
two equations are normalized, respectively, by 



^ P(AAge) d(AAge) = 1, 
and 

^ P(AZ)d(AZ) = 1. 



(5) 



(6) 



To see the fits of the two distributions clearly, we plot the distri- 
butions derived from fitting {B - V) and {B — K) colors and the 
distributions calculated by Eqs. (3) and (4) in Fig. 12. In prac- 
tice, these distributions can possibly be used to coiTect for the 
deviations of stellar ages and metallicities determined by colors 
and get some reliable estimations for the two parameters. 
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Fig. 11. Distributions of the differences between stellar parame- 
ters determined by fitting color pairs and those of the DSPs of 
CSPs (AAge = ff - fi, AZ = Zf - Zi). 
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Fig. 12. The fitted distributions of AAge and AZ. AAge and AZ 
have the same meanings as in Fig. 11. Solid lines show the dis- 
tributions derived from a sample of CSPs and dashed lines show 
the distributions calculated by Eqs. (3) and (4). Panels (a) and (b) 
show the distributions for, respectively, the differences in stellar 
age and in metalUcity determined by [(B -V) ,{B- K)]. 
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5.2. Correcting the deviations of stellar ages and metallicities 

We tried to correct for the deviations in stellar ages and metal- 
licities determined by fitting a pair of colors for a sample of 2 
000 CSPs with f 1 > 3 Gyr The sample is generated by the same 
method as in Sect. 4 and all of their UBVRIJHK colors are cal- 
culated. We get the best-fitted stellar ages and metallicities (ff 
and Zf) for CSPs by fitting their (B - V) and (B - K) colors. 
The distributions of the differences between the best-fitted pa- 
rameters and those of their DSPs (AAge = tf - t\ and A(Fe/H) 
= [Z/H]f - [Z/H]i) are then calculated, and are shown in panels 
(a) and (b) of Fig. 13. We see that some stellar ages and metal- 
licities younger and richer than those of the DSPs of CSPs are 
determined by color pairs. This suggests that without any cor- 
rections, it would be difficult to get further reliable results based 
on the stellar ages and metallicities determined by color pairs. 

To correct for the deviations in stellar ages and metallici- 
ties derived from colors, we tried two methods. In the first, we 
corrected the stellar ages and metallicities derived from colors 
using two maximum-likelihood values, i.e., -0.8 Gyr and 0.0002 
for stellar age and metallicity, respectively. The results show that 
about 50% CSPs have stellar ages and metallicities close to those 
of the DSPs of these populations after the correction (AAge < 1 
Gyr, A([Fe/H]) < 0.05). In the second, we corrected stellar ages 
with the distribution of AAge, which is shown in Fig. 12. Each 
corrected age was obtained by subtracting a random value that 
fits to the distribution shown in panel (a) of Fig. 12 from the best- 
fitted age, ff . Stellar metallicities are corrected in a similar way, 
using the distribution shown in panel (b) of Fig. 12. We call the 
correction "distribution correction". For convenience, we use the 
distributions obtained by this work rather than ones calculated 
by Eqs. (3) and (4). The distribution of the differences between 
corrected stellar ages and those of their corresponding DSPs can 
be seen in the left column of Fig. 18, and Fig. 12 is useful for 
comparison. We see clearly that the stellar ages corrected by us 
distribute homogeneously around those of their DSPs. A similar 
correction is also used to correct the stellar metallicities and the 
results shown in the right column of Fig. 18. As can be seen, 
the metallicities of CSPs are homogeneously distributed around 
those of their DSPs after the correction (see Fig. 12 for com- 
parison). Therefore, the distribution correction can help us to 
get stellar ages and metallicities distributed around those of the 
DSPs of CSPs. To see the effects of the correction more clearly, 
we plot the results in another way (Fig. 13). In the plot, uncor- 
rected stellar ages, log(ff/yr), uncorrected stellar metallicities, 
[Z/H]f=log(Zf/Zo), corrected stellar ages, log(fc/yr), and cor- 
rected stellar metallicities, [Z/H]c=log(Zc/Zo) of 2 000 CSPs 
are plotted versus those of their DSPs (log(ri/yr) and [Z/H]i). 
Most stellar ages and metallicities corrected by us fit to those of 
their DSPs when taking possible uncertainties in stellar age (1.5 
Gyr) and metallicity (0.005) into account. Therefore, the aver- 
ages of stellar ages and metallicities of a sample of galaxies can 
be estimated in this way, as can their distributions of stellar ages 
and metallicities. 

We apply the correction to the results of 53 elliptical galax- 
ies, whose stellar ages and metallicities are estimated by {B - V) 
and {B - K) colors by Li et al. (2007). The distributions of un- 
corrected and corrected stellar parameters are plotted in Fig. 14. 
After the correction, stellar populations of these galaxies, which 
average about 6 Gyr and 0.0375 in stellar age and metallicity, 
respectively, come to be older and less metal rich than the uncor- 
rected ones. The results seem to be closer to the typical values 
of early type galaxies that are determined by Lick/IDS indices 



or spectra (see, e.g., Gallazzi et al. I2005I I. although these popu- 
lations are somewhat younger and more metal rich. 

To check the reliability of stellar ages and metallicities 
measured by colors, we fit the two parameters for 61 early 
type galaxies respectively by two line-strength indices (HyS and 
[MgFe]= VMgb X (Fe5270 + Fe5335)/2 , see Gonzalez [1993] ) 
and two colors {B -V and B - K), and compare the two kinds 
of results. Although there is a new definition of [MgFe], in 
which Fe5270 and Fe5335 are weighted as 0.72 and 0.28, re- 
spectively (Thomas et al.'2003), we take the the above definition, 
because we only have the data for < Fe > = (Fe5270+Fe5335)/2. 
The Lick indices of the sample galaxies are selected from 
Sil'chenko ( |2006] l, Trager ('20U0'), and Thom as et a l. (,2005i ), 
which can also be found from Kuntschner et al. (I2006I I. The pho- 
tometry data and redshifts of these galaxies are obtained from 
the NASA/IPAC Extragalactic Database (NED) and HyperLeda 
database (http://leda.univ-lyonl.fr/), respectively. The K-bwd 
magnitudes are calculated from /Tj-band data and take the k- 
correction and galactic extinction into account. Here, we take 
the methods presented by Bessell (2005 ) & Girardi et al. (120031 ) 
for transforming K^-bwd data into K-bwd and ^-correcting the 
data. The galactic extinction is calculated using the calculator 
supplied by NED. 

We list the stellar ages and metallicities fitted by the two 
methods in Table 3. Note that we only show 18 galaxies that have 
best-fitted stellar ages (fuck and fcoior) smaller than 15 Gyr here, 
according to the age of the universe (see, e.g., Shafieloo et al. 
2006 ). The results in the table were not corrected for the effects 
of YSPs, so that we can easily estimate the difference between 
the stellar-population parameters determined by colors and those 
by line-strength indices. We did not correct for the light travel 
time either, because the correction is not important for nearby 
galaxies. In Fig. 15, we compare the stellar parameters deter- 
mined by colors and corrected via "distribution correction" with 
those determined by line strength indices. Note that the stellar 
ages and metallicities derived from colors are corrected in the 
whole sample of 61 galaxies. Possible uncertainties (1.5 Gyr in 
age and 0.005 in metallicity) are shown in Fig. 15. We see that 
the two kinds of results are consistent, taking uncertainties about 
2.5 Gyr in age and 0.01 in metallicity. 

When we try to compare the stellar ages and metallicities 
estimated respectively by a color pair, {{B-V), (B-K)}, and four 
independent colors, (B - V), {B - K), (/ - H), and {B - I), we 
obtain similar results. The colors of galaxies are selected from 
Michard (120051 ) . The best-fitted stellar ages and metallicities are 
found by minimizing the function 



(7) 



where C,- is the /th observed color, C. the color corresponding 
to C, for the jth SSP, and n (2 or 4) the number of colors used for 
fitting. In the fitting, we did not take the observational uncertain- 
ties into account, because there is no reliable data. The possible 
ranges of stellar ages and metallicities determined by a color pair 
and by four independent colors are given at 1 cr (68.3%) and 2 
cr (95.4%) confidence levels. The confidence levels can be es- 
timated via a chi-square test, taking a typical uncertainty (0.02 
mag) for each color. The detailed comparison are seen in Fig. 16, 
in which the results of four elliptical galaxies are shown. From 
the figure, we see that the possible stellar ages and metallicities 
are grouped along the line of the well-known age-metallicity 
degeneracy (Worthev ll994| ), and the possible ranges determined 
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Fig. 13. Plots of uncorrected and corrected stellar ages and 
metallicities versus those of their corresponding DSPs. Panels 
(a) and (b) show the plot of the uncorrected stellar parameters 
(log(ff/yr) and [Z/H]f) and those of their corresponding DSPs 
(log(fi/yr) and [Z/H]i). The uncorrected stellar ages and metal- 
licities are determined by {B - V) and (B - K) colors. Panels (c) 
and (d) show the similar relations for distribution-corrected stel- 
lar ages and metallicities (log(f(;/yr) and [Z/HJc). Dashed lines 
in panels (a) and (c) show a + 1.5 Gyr spread about the unity 
(solid) line for stellar ages. In panels (b) and (d), they show and 
a + 0.005 spread about the unity (soUd) hne for stellar metalhc- 
ities. 



by the two methods are similar. This can be interpreted as there 
being no color sensitive only to DSP or YSP among the four col- 
ors. One can refer to Li et al. (I2007I I. in which a principal com- 
ponent analysis (PCA) of colors for SSPs with various ages was 
shown. However, most of the best-fitted stellar ages and metal- 
licities determined by four colors are respectively younger and 
more metal rich than those determined by a color pair, as we see. 
This suggests that the YSPs in galaxies affect the determination 
of DSPs of galaxies more strongly when we measure stellar ages 
and metalhcities by fitting a few independent colors. 

5.3. Effect of the calculation of mass fractions of young 
populations 

We took 0.5 and 3.02 for Fq and t, respectively, to calculate the 
mass fractions of YSPs of CSPs via Eq. (2) in our standard case, 
but different values may lead to different distributions of devia- 
tions in stellar ages and metallicities. Here we test how Fq and 
T change the distributions of deviations in the stelllar-population 
parameters estimated by two colors. Given t = 3.02, we test the 
effects of F() via four values, i.e., 0.2, 0.3, 0.4, and 0.5. The dis- 
tributions relating to different values of Fq are plotted in panels 
(a) and (b) of Fig. 17. Similarly, given Fo - 0.5, we test the ef- 
fects of T, taking 2.0, 3.02, 5.0, and 10.0. The results are shown 
in panels (c) and (d) of Fig. 17. In addition, we test the distribu- 
tions by taking fixed values (5%, 10%, and 15%) for F2, which 
are shown in panels (e) and (f) of Fig. 17. As we see, high val- 
ues of Fq or T will lead to more populations with big deviations 
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Fig. 14. Distributions of uncorrected and corrected stellar ages 
and metallicities of 53 elliptical galaxies. Panels (a) and (b) 
show the distributions of the uncorrected results (log(ff/yr) and 
[Z/H]f) while panels (c) and (d) show those for the corrected 
results (log(fc/yr) and [Z/H]c). 

Table 3. Stellar ages and metallicities fitted respectively by line 
strength indices (H/? and [MgFe]) and colors (B -V and B - K). 



Galaxy 


log(fiick/yr) 


[Z/H]iick 


log(fcolor/yr) 


[Z/H]color 


NGC0080 


10.0019 


0.1072 


9.6750 


0.2707 


NGC2681 


9.5119 


-0.4815 


9.3900 


-0.1707 


NGC2950 


9.5119 


0.3627 


9.4950 


0.1917 


NGC3166 


9.5119 


-0.0066 


9.7650 


-0.0809 


NGC3377 


9.5969 


0.1553 


9.9450 


-0.4377 


NGC3384 


9.9469 


-0.0022 


9.7600 


0.0881 


NGC3941 


10.1269 


-0.1612 


9.7700 


0.0374 


NGC4026 


9.9769 


0.0000 


9.6600 


0.1492 


NGC4150 


9.7769 


-0.5302 


9.8550 


-0.5768 


NGC4179 


9.9169 


0.1847 


9.5500 


0.2055 


NGC4350 


9.9769 


0.1173 


9.6700 


0.2455 


NGC4570 


9.9519 


0.1945 


9.6600 


0.2201 


NGC5308 


9.9369 


0.2822 


9.5300 


0.2240 


NGC5574 


9.5119 


-0.0757 


9.5500 


-0.1397 


NGC7013 


9.7319 


-0.1457 


9.4300 


-0.6289 


NGC7280 


9.9269 


-0.2147 


9.4700 


0.1055 


NGC7457 


9.6519 


-0.1457 


9.9350 


-0.4685 


NGC7743 


9.6869 


-0.1367 


9.4750 


0.0128 



Note : Suffixes "lick" and "color" represent the results fitted by line 
strength indices and those by colors, respectively. 



in two stellar-population parameters. However, the distributions 
of deviations in the two stellar-population parameters are similar 
for different formulations of F2. We tried to correct for the de- 
viations using the distributions obtained in the standard case, by 
a few samples including 2 000 CSPs generated by different for- 
mulations of F2. The distributions of differences between two 
corrected stellar-population parameters and those of the DSPs 
are shown in Fig. 18. In the test, samples of CSPs are generated 
by taking different formulations. As we see, the stellar ages and 
metallicities distribute around those of DSPs, after the distribu- 
tion correction. Therefore, even if the star formation histories are 
different, the averages of the ages and metallicities of the DSPs 
of galaxies can be estimated from their colors, by using the the 
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Fig. 16. Comparison of the stellar parameters (age and metallicity) determined by a color pair (B-V and B-K) with those determined 
by four independent colors {B - V, B - K, B - I, I - H). Squares and triangles show the best-fit stellar ages and metallicities that 
are determined by a color pair and four independent colors, respectively. Gray and dark lines show the 1 cr (solid) and 2 cr (dashed) 
ranges of stellar ages and metallicities for the two-color method and multi-color method, respectively. 




0.05 
A([Fe/H]) 

Fig. 17. Distributions of the deviations in stellar ages and metallicities determined by color pairs for different calculations of the 
mass fractions of YSPs of composite star systems. Panels (a, b), (c, d), and (e, f) show for varying Fq, t, F2, respectively. In panels 
(a), (b), (c), and (d), F2 are calculated by Eq. (2). Our standard case takes Fo=0.5 and t=3.02 for Eq. (2). 



distributions obtained in the standard case to correct the results, 
as can the distributions of the two stellar-population parameters. 



6. Conclusions and discussions 

We first investigated the effects on colors of a composite stellar 
population (CSP) of three parameters (age, t2, mass fraction, F2, 
and metallicity, [Z/H]2) of the young stellar population (YSP). 
Then we studied the potential distributions of the deviations of 



Zhongmu Li and Zhanwen Han: Effects of YSPs on determination of stellar populations 



9 



CM 

d 



c 
o 



2 6 



d 



o 



CM 

d 

c 
o 

.2d 



[i hiT^fpf 1 1 1 


. (a): 
iiVi'frT'l-i- j I- 


1 1 ' 

> i .1 .j^.-T^ 


1 ' 1 ' . 
t] 

I^V-^-t^i — '■ 


- 


; (c): 

Ji : 

• J 

1 1 1 Tyt^i^ 't I- 


r 

; ■ 1 ■ ■ ■jt.f-r^ 




Li ^^^^T|•T^"l 1 1 

_ 1 


(e): 

ll ■ 

1 _ 


r 

L , J 


(f): 



-10 - 



5 5 
Mge (Gyr) 



10 



-0.2 0.2 
A([Fe/H]) 



CM 

d 



CM 

d 



CM 

d 



Fig. 18. Distributions of the deviations in corrected stellar ages and metallicities from those of their coiTesponding DSPs. Solid lines 
show for 2 000 CSPs with F2 calculated by Eq. (2) with Fo = 0.5 and t - 3.02. The other samples of CSPs in Panels (a) and (b) are 
generated by taking Fo - 0.2 (dash-dot-dot-dot), 0.3 (gray solid), 0.4 (dotted), and r = 3.02 for Eq. (2). The ones in Panels (c) and 
(d) are generated by taking t - 2 (dash-dot-dot-dot), 5 (gray solid), 10 (dotted) and Fq - 0.5 for Eq. (2). Those in Panels (e) and (f) 
are generated by taking F2 = 5% (dash-dot-dot-dot), 10% (gray solid), 15% (dotted), respectively. 
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Fig. 15. Stellar ages and metaUicities fitted by [(B - V), (B - K)] 
versus those fitted by Hfi and [MgFe]. Dashed lines show a + 2.5 
Gyr spread about the unity (solid) line for stellar ages in panel 
(a), while they show a + 0.01 spread in stellar metallicities in 
panel (b). Error bars show uncertainties of 1.5 Gyr and 0.005 in 
stellar age and metalhcity, respectively. 



stellar ages and metallicities determined by color pairs, based on 
a set of CSPs that are rebuilt via SSPs of the BC03. 

Our results show that the age (fa) and mass fraction (F2) 
of the YSP affect colors of the system significantly, while the 
metallicity of YSP (Z2 or [Z/H]2) affect the system's colors more 
weakly. Thus if there were YSPs in galaxies and they were not 
too old (e.g., f2 < 5 Gyr), the colors of galaxies should be sig- 
nificantly different from those of their DSPs. This means that it 
is difficult to get the accurate ages and metallicities of the DSPs 
of galaxies by comparing colors of galaxies to those of SSPs di- 
rectly. In addition, the effects of the age (f2) and mass fraction 
(F2) of the YSP are found to be degenerate: The smaller the age 
or greater mass fraction, the bigger the deviation in the colors 
of a CSP compared to those of the DSP of the system. But, as a 
whole, the effects of t2 on the colors of star systems are stronger 



than those of F2. Our results also show that the effects of the 
metallicity of the YSP (Z2 or [Z/H]2) on the colors of a system 
are much weaker than those of f2 and F2- 

Furthermore, we give possible distributions for the devia- 
tions in the stellar ages and metallicities determined by color 
pairs in this work, which can help us to get reliable estimations 
of the averages and distributions of the stellar ages and metallic- 
ities of a sample of galaxies from their colors. This is possibly 
the most important result of this work. According to our results, 
the stellar ages and metallicities measured by color pairs, e.g., 
[(B - V), (V - K)], [(B - V), (B - K)l and [(B - V), (I - H)l 
are about 2.14 Gyr younger, while 0.0027 more metal rich than 
those of the DSPs of CSPs on average. The differences between 
the colors of CSPs and those of their DSPs were also studied in 
this work. 

When we tried to study the stellar ages and metallicities de- 
termined respectively by a color pair and four independent col- 
ors, we got similar ranges for the two stellar-population parame- 
ters. The results also suggest that YSPs affect the determination 
of the ages and metallicities of DSPs more strongly when we 
estimate the two stellar-population parameters by a few inde- 
pendent colors rather than a color pair. 

Our results can be used widely, although we took the BC03 
models for the study. The reason is that various models usually 
predict similar colors for the same stellar population. For ex- 
ample, (V - K) colors predicted by PEGASE models are only 
significantly different from those predicted by the BC03 models 
when the stellar age is between 0.06 and 0.32 Gyr, but (B - V) 
colors predicted by the two models are very similar for all SSPs 
(Bruzual & Chariot |2003]l. 

In addition, the results can help us not only estimate the 
stellar ages and metallicities but also study the fundamental 
relations, such as the fundamental plane (FP) (Djorgovski & 
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Davis 119871 ) and Kormendy relation (Kormendy 1 19771 1 of el- 
liptical galaxies and the bulges of lenticulars, which are power- 
ful distance indicators for galaxies and important constraints on 
galaxy formation and evolution. According to our results, YSPs 
increase B luminosity, and, correspondingly, B surface bright- 
ness, compared to that of the DSPs of galaxies. In other words, 
YSPs affect the light of galaxies significantly. However, YSPs 
are unimportant for affecting the dynamic parameters of galax- 
ies, as they only contain a small fraction of the stellar mass. 
Therefore, it is necessary to take the effects of the YSPs of 
galaxies into account when we study the fundamental relations 
relating to the dynamic parameters of galaxies, e.g., the rela- 
tions transformed from the virial theorem. The FP of early type 
galaxies is just one of these kinds of relations. Using the stellar 
populations of 18 galaxies, which are measured from two colors 
(B-V and B-K) and two line strength indices (HyS and [MgFe]), 
we calculated the difference between B surface brightness of a 
CSP and that of its DSP. The results show that the difference 
between B surface brightness of a CSP and its DSP can be as 
large as 1 mag, with an average of 0.33 mag. Therefore, taking 
only the DSPs in early type galaxies into account, the FP of early 
type galaxies will change: A pair of effective radius ( rg) and ve- 
locity dispersion (cr) will correspond to fainter B surface bright- 
ness compared to the FP obtained before. In this case, we will 
measure distance moduli smaller by about 0.33 mag for galaxies 
than those determined via the FP obtained before. Furthermore, 
if we correct for the effects of YSPs, the FP obtained in vari- 
ous bands may change. Possibly, the differences among FPs ob- 
tained in various bands may decrease, compared to the results of 
Bernardi et al. (2003 ). The Kormendy relation (Kormendv l 1 9771 l 
will change and affect the determination of the distances of early 
type galaxies with a trend similar to the FP if we take the effects 
of YSPs into account. A given effective radius corresponds to 
fainter B surface brightness and then give us a smaller distance 
than the previous relation. Moreover, we find that the average of 
the difference between B surface brightness of a CSP and that of 
its DSP, i.e., 0.33 mag, is close to the typical scatter of the FP 
and Kormendy relation of elliptical galaxies and the bulges of 
lenticulars (0.3-0.5 mag) (see, e.g., Falcon-Barroso et al. 20021 
Djorgovski & Davis flMTt Kormendv 119771 Reda et al. I2005] ). 
It implies that YSPs in early type galaxies may contribute to the 
scatter of the two fundamental relations. Thus the scatter in the 
two fundamental relations may be smaller if we take the effects 
of YSPs into account. 

In conclusion, the results of the work can help us to get stel- 
lar ages and metallicities close to those of the DSPs of galaxies 
via photometry data and to understand the fundamental relations 
of early type galaxies better. The main results of this work are 
available in electronic form at the CDS. 
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